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ABSTRACT 

We use the newest generation of the Starburst99/Mappings code to generate an extensive suite of 
models to facilitate detailed studies of star- forming galaxies and their ISM properties, particularly 
at low metallicities. The new models used include a rigorous treatment of metal opacities in the 
population synthesis modeling and more detailed dust physics in the photoionization code. These 
models span a wide range of physical parameters including metallicity, ionization parameter, and 
the adoption of both a instantaneous burst and continuous star formation history. We examine the 
agreement between our models and local {z < 0.1) star- forming galaxy populations from several large 
datasets, including the Sloan Digital Sky Survey, the Nearby Field Galaxy Survey, and samples of blue 
compact galaxies and metal-poor galaxies. We find that models adopting a continuous star formation 
history reproduce the metallicity-sensitive line ratios observed in the local population of star-forming 
galaxies, including the low-metallicity sample. However, we find that the current codes generate an 
insufficiently hard ionizing radiation field, leading to deficiencies in the [SII] fiuxes produced by the 
models. We consider the advantages and short-comings of this suite of models, and discuss future 
work and improvements that can be applied to the modeling of star- forming galaxies. 
Subject headings: 



1. INTRODUCTION 

Robust analysis of star-forming galaxy emission-line 
spectra can provide constraints on key physical param- 
eters of the ionizing radiation field and the interstellar 
medium (ISM). The star formation rates (SFRs) of these 
galaxies can be estimated from luminosities of the Ha 
line (Hunter & Gallagher 1986, Kennicutt 1998, Bicker 
& Fritze-v. Alvensleben 2005, Kewley et al. 2007) or 
the [Oil] line (Gallagher et al. 1989, Kennicutt 1998, 
Rosa-Gonzalez et al. 2002, Kewley et al. 2004, Mous- 
takas et al. 2006), and multiple studies have examined 
the use of the H/3 equivalent width to estimate the age 
of the young stellar population in galaxies (Schaerer & 
Vacca 1998, Gonzalez Delgado et al. 1999, Fernandes 
et al. 2003, Martin-Manjon et al. 2008). Metallicity is 
another critical parameter in the study of star-forming 
galaxy ISM environments, shedding light on star forma- 
tion histories and subsequent chemical evolution. A va- 
riety of optical emission line ratio diagnostics have been 
presented and employed to determine metallicities, in- 
cluding measuring the [OIII] line ratios to determine elec- 
tron temperatures, and therefore abundances (e.g. Peim- 
bert 1967, Garnett 1992), and use of the [NII]A6584/Ha 
and [Nil] A6584/ [OIII] A5007 ratio diagnostics from Pet- 
tini k Pagel (2004), the [NII]A6584/[OII]A3727 ratio 
from Kewley & Dopita (2002), and the ([OIII]A5007 -I- 
[OIII]A4959 + [OII]A3727)/H/3 (R23) ratio (e.g., Pagel et 
al. 1979, McGaugh 1991, Zaritsky et al. 1994, Pilyugin 
2000, Chariot & Longhctti 2001, Kewley & Dopita 2002, 
Kobulnicky & Kewley 2004). 

Baldwin et al. (1981) present the technique of plot- 
ting optical emission line ratios, such as [NII]A6584/Hq; 
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VS. [OIII]A5007/H/?, on a series of diagnostic diagrams 
to separate extragalactic objects according to their pri- 
mary excitation mechanisms. These and other diagrams 
were later used by Veilleux & Osterbrock (1987) to de- 
rive a semi-empirical classification scheme to distinguish 
between star-forming galaxies and active galactic nuclei. 
Kewley et al. (2001) used these same diagrams to derive 
a purely theoretically classification scheme, which was 
later extended by Kauffmann et al. (2003) and Stasin- 
ska et al. (2006) using data from the Sloan Digital Sky 
Survey. Emission line diagnostic diagrams can also be 
used to probe the shape of the far-ultraviolet ionizing 
spectrum of a galaxy (Dopita et al. 2000, Kewley et al. 
2001, Kewley & Dopita 2002). As a result, diagnostic 
ratio grids are often employed to test the agreement of 
stellar population synthesis and photoionization model 
grids with emission-line ratios measured in observations 
of star-forming galaxies. These comparisons effectively 
illustrate the evolution and improvements of such mod- 
els, and are also useful in highlighting the shortcomings 
of different grids and the challenges faced in modeling 
these galaxies and their emission spectra. 

Kewley et al. (2001) use both the Starburst99 (Lei- 
therer et al. 1999) and Pegase (Fioc & Rocca-Volmerage 
1997) evolutionary synthesis models in conjunction with 
the Mappings III photoionization code (Binette et al. 
1985, Sutherland & Dopita 1993) to compute photoion- 
ization models that they compare to a sample of 157 
warm IR starburst galaxies on a variety of optical emis- 
sion line ratio diagnostic diagrams. They find that as- 
suming a continuous star formation history is more realis- 
tic than the assumption of a single instantaneous burst. 
Moy et al. (2001) find a similar result; using the Pe- 
gase evolutionary synthesis models and the CLOUDY 
photoionization code (Ferland 1996) they model young 
stellar populations along with underlying stellar popula- 
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tions, and find that underlying popiilations with contin- 
uous star formation histories are more compatible with 
observed emission line ratios in starburst galaxies. Fer- 
nandes et al. (2003) also find support for modeling un- 
derlying stellar populations in their population synthesis 
analysis of starburst and HII galaxies, using the Bruzual 
& Chariot (1993) evolutionary synthesis code. 

Kewley et al. (2001) do, however, find that the ionizing 
spectra produced by the Starburst99 models are not hard 
enough in the far ultraviolet (FUV) region of the spec- 
trum to reproduce the observed line ratios, and propose 
that including the effects of continuum metal opacities 
in stellar atmospheres should be a way of improving the 
models' accuracy. This deficiency in the FUV ionizing 
spectrum is noted in Panuzzo et al. (2003) and Magris 
et al. (2003). StarburstQQ models at the time used older 
stellar atmosphere models that do not include treatments 
of opacities beyond hydrogen and helium (Schmutz et al. 
1992), or else are not complete at the high temperatures 
or masses that prove critical when modeling the FUV 
ionizing radiation field (Rauch 1997, Lejeune et al. 1997, 
1998). 

Another challenge that is common among previous 
stellar population synthesis and photoionization models 
is a difficulty in modeling the environments of metal-poor 
galaxies (e.g. Fernandes et al. 2003, Dopita et al. 2006, 
Martin-Manjon et al. 2008). Dopita et al. (2006) use the 
latest Starburst99 code from Vazquez & Leitherer (2005) 
and the Mappings III photoionization code to generate 
models of isobaric dusty HII regions, as well as integrated 
galaxy spectra that model the galaxies as composites 
of multiple HII regions. Dopita et al. (2006) removes 
the ionization parameter as a free parameter, fixing the 
ISM pressure and central cluster mass and instead al- 
lowing ionization parameter to vary with time, pointing 
out that ionization parameter also varies with metallic- 
ity. Dopita et al. (2006) find that these models still do 
not reproduce the observed emission line diagnostic ra- 
tios of lower-metallicity galaxies {Z < OAZq). They 
also note that the ionizing fluxes are not hard enough 
to agree with the integrated spectra. Martin-Manjon et 
al. (2008, 2009) take a different approach, using sets of 
stellar yields from Gavilan et al. (2005) to model the 
chemical evolution of HII galaxies, allowing metallicity 
to evolve with age rather than generating a grid with 
a range of fixed metallicities. They also use the newer 
generation of CLOUDY (Ferland et al. 1998) to model 
environments that have undergone multiple bursts of star 
formation (Martin-Manjon et al. 2009), but still find that 
these models don't account for the most metal-deficient 
HII galaxies in their sample. 

Metal-poor galaxies are an important avenue of study 
in their own right. With a gas-phase oxygen abundance 
of log(0/H) + 12 < 7.9 (Yin et al. 2007), they offer 
a relatively pristine pre-enchriment ISM in which star- 
formation, as well as current and previous episodes of 
enrichment, can be examined (Brown et al. 2008). These 
galaxies also pose challenges to scenarios seeking to re- 
produce their metal-poor environments through a variety 
of evolutionary mechanisms (see Kewley et al. 2007 and 
references therein). It has also recently been proposed 
that the host galaxies of long-duration gamma-ray bursts 
belong to the metal-poor galaxy population (Stanek et 
al. 2006, Fruchter et al. 2006, Kewley et al. 2007, Modjaz 



et al. 2008). Proper modeling of these galaxies is criti- 
cal to our understanding of the star formation processes 
and evolving stellar populations present in such environ- 
ments. 

Here we present models generated using the Vazquez & 
Leitherer (2005) Starburst99 stellar population synthesis 
code and the latest generation of the Mappings III code, 
with recent improvements that include a more rigorous 
treatment of dust (Groves et al. 2004). Our models are 
tailored towards addressing the difficulties found in past 
work with producing a harder FUV spectrum; we do this 
by adopting the WMBASIC stellar atmosphere models 
of Pauldrach et al. (2001) and the CMFGEN HiUier & 
Miller (1998) atmospheres, both of which include the rig- 
orous treatments of continuum metal opacities that were 
suggested in Kewley et al. (2001). Our grid was also de- 
signed to precisely model the emission line flux of the 
local galaxy population, including low-metallicity galax- 
ies, adopting the full range of metallicities available in 
the evolutionary tracks published by the Geneva group 
(Schaller et al. 1992; Schaerer et al. 1993a, 1993b; Char- 
bonnel et al. 1993). 

In this paper we describe our new grid of stellar pop- 
ulation synthesis models, detailing the inputs and free 
parameters that we adopt when generating the grid and 
examining the ionizing FUV spectra that arc produced 
by Starburst99 in detail (§ 2). We consider how these 
spectra affect the behavior of a variety of optical emis- 
sion line diagnostics (§ 3). With these diagnostic ratios, 
we generate a series of optical emission line diagnostic 
diagrams, and compare our model grids to spectra of a 
variety of nearby (z < 0.1) galaxy populations to assess 
their agreement (§ 4) . Finally, we consider the results of 
these comparisons and discuss potential future work in 
this area (§5). 

2. STARBURST99/MAPPINGS III MODEL GRIDS 
2.1. Model Grid Parameters 

To model our sample of galaxies we have used the Star- 

burst99 code (Leitherer et al. 1999, Vazquez & Leitherer 
2005) to generate theoretical spectral energy distribu- 
tions (SEDs), which in turn were used in the Mappings 
HI photoionization models to produce model galaxy 
spectra that could be compared to our observations. 

Starburst99 is an evolutionary synthesis code that can 
be used to generate synthetic ionizing far-ultraviolet 
(FUV) radiation spectra as a fimction of metallicity, star 
formation history, and the age and evolution of the stel- 
lar populations. These populations are produced by use 
of model stellar atmospheres and spectra along with evo- 
lutionary tracks for massive stars. For this work, we have 
used a Salpeter initial mass fimction (Salpetcr 1955) with 
a IOOMq upper mass boundary, along with Starburst99's 
Pauldrach/Hillicr model atmospheres. These employ the 
WMBASIC wind models of Pauldrach et al. (200i) for 
younger ages when O stars dominate the luminosity (< 
3 Myr), and the CMFGEN HiUier & Miller (1998) atmo- 
spheres for later ages at which Wolf-Rayet stars are domi- 
nant. This differs from the Starburst99 models presented 
in Dopita et al. (2000) and Kewley et al. (2001), which 
adopt the plane-parallel Lejemie et al. (1997) grid of at- 
mosphere models along with the Schmutz et al. (1992) 
extended model atmospheres for stars with higher winds. 
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The Schmutz et al. (1992) models, which include the crit- 
ical Wolf-Rayct phase, only include continuous opacities 
for hydrogen and helium, neglecting what are expected 
to be considerable effects from continuum metal opac- 
ities. Kewley et al. (2001) suggest this as a potential 
shortcoming in their models, proposing that the inclu- 
sion of continuum metal opacities will result in a fraction 
of higher-energy radiation being absorbed and reemit- 
ted at lower energies, in the region of the FUV spec- 
trum that is responsibly for ionizing the optical emission 
lines. The Pauldrach et al. (2001) and Hillier & Miller 
(1998) model atmospheres address this shortcoming by 
including rigorous non-LTE treatments of metal opaci- 
ties. In conjunction with the Pauldrach/Hillicr atmo- 
spheres we also adopt two different sets of evolutionary 
tracks produced by the Geneva group, and consider the 
particulars of these tracks' mass loss rates (see Section 
§ 2.2). Starburst99 generates the final synthetic FUV 
spectrum output through use of the isochronc synthesis 
method first introduced by Chariot & Bruzual (1991), fit- 
ting isochrones to the evolutionary tracks across different 
masses rather than discretely assigning stellar mass bins 
to specific tracks. 

The resulting FUV spectrum is then taken as input by 
the Mappings III code. The Mappings shock and pho- 
toionzation code was originally developed by Bincttc et 
al. (1985), improved by Sutherland & Dopita (1993), and 
used in Dopita et al. (2000) and Kewley et al. (2001). 
The recent improvements to the Mappings code, used in 
Dopita et al. (2006) and Snijders et al. (2007), include a 
more sophisticated treatment of dust (Groves et al. 2004) 
that is adopted in our models, which properly treats the 
effects of absorption, charging, and photoelectric heating 
by the grains - for a more detailed discussion, see Groves 
et al. (2004) and Snijders et al. (2007). Taking the syn- 
thetic ionizing FUV spectrum output of Starburst99 and 
an adopted nebular geometry model, which we assume 
to be plane-parallel, we select a variety of electron den- 
sities and ionization parameters as inputs. Using these 
parameters in Mappings III we computed a complete grid 
of plane-parallel isobaric photoionization models. 

When generating our Starburst99/Mappings III model 
spectra, we adopted a broad grid of input parameters to 
facilitate comparison with a wide range of galaxy sam- 
ples: 

Star Formation History (SFH): We model both a zero- 
age instantaneous burst of star formation, with a fixed 
mass of lO^Af0, and a continuous SFH where the star 
formation rate (SFR) is constant at a rate of 1 Mq per 
year, starting from an initial time and assuming a stellar 
population that is large enough to fully sample the IMF 
at high masses. 

Metallicity: We adopt the full range of metallicities 
available from the evolutionary tracks of the Geneva 
group, which includes five metallicities of z = 0.001 
{Z = 0.05^0), z = 0.004 {Z = 0.2Zq), z = 0.008 
{Z = 0.4Zq), z = 0.02 {Z = Zq), and z = 0.04 
(Z = 2Zq). 

Evolutionary Tracks: We adopt the two evolutionary 
tracks of the Geneva group that are currently available 
in Starburst99: the Geneva "Standard" mass loss tracks, 
and the Geneva "High" mass loss tracks. The differences 
in these tracks are discussed in more detail in Section 
§ 2.2. 



Age: We generate models ranging in age from to 5 
Myr in the case of an instantaneous burst star formation 
history in 0.5 Myr increments. For the continuous SFH 
models we adopt a constant age of 5 Myr, the age at 
which a continuous SFH stellar population reaches equi- 
librium (Kewley et al. 2001). Further discussion of the 
age range of these model grids can be found in § 2.3. 

Ionization parameter. The ionization parameter q (cm 
s~"^) can be thought of as the maximum velocity possible 
for an ionization front being driven by the local radiation 
field. The value itself is calculated for the inner surface of 
the nebula. By this definition q relates to the dimension- 
less ionization parameter lA hy U = q/ c. Rigby & Rieke 
(2004) find a range of —3 < log U < —1.5 for the dimen- 
sionless ionization parameter in local starburst galaxies. 
In our model grid, we adopted seven different values for 
our ionization parameter q, where q = Ix lO"^. 2 x 10^, 4 x 
10^, 8 X 10^, 1 X 10^, 2 X 10*, and 4 x 10* cm s^^). These 
values correspond to dimensionless ionization parameters 
of log U w -3.5, -3.2, -2.9, -2.6, -2.5, -2.2, and -1.9, 
respectively. While these are slightly lower than the log 
U values found in Rigby & Rieke (2004), they are similar 
to the range of ionization parameters adopted in Kewley 
(2001) and Snijders et al. (2007). 

Electron density: Dopita et al. (2000) find — 10 
cm~^ to be typical of giant extragalactic HII regions, 
while Dopita et al. (2006) constrain the electron density 
in their models of HII regions to rie ^ 100 cm~^. In 
addition, Kewley et al. (2001) find an average electron 
density rig = 350 cm~'^ for the Kewley et al. (2000) sam- 
ple of warm infrared starburst galaxies. Since we wish to 
compare our models to normal and low-metallicity star- 
forming galaxies, which are expected to have lower rie 
than those found in the gas-rich warm infrared galaxies, 
we adopt two different electron densities in our full model 
grid, adopting both rie = 10 cm~^ and Ue = 100 cm~^. 
We assume an isobaric density structure for these mod- 
els, and thus is specified by the dimensionless pres- 
sure/mean temperature ratio. We found that the lower 
Ue = 10 cm~^ produced only a very slight decrease in 
flux for the more g-sensitive lines at higher metallicities. 
For the remainder of this paper, we present results that 
adopt Ue = 100 cm~^, following the findings of Kewley 
et al. (2001) for starburst galaxies. 

2.2. Stellar Evolutionary Tracks 

Currently, Starburst99 includes two different evolu- 
tionary tracks produced by the Geneva group, which dif- 
fer primarily in their treatment of mass loss rates for 
massive stars. Mass loss rates are a critical parameter 
when considering the contributions of massive stars to 
ISM enrichment (Maeder & Conti 1994). 

The "standard" (STD) mass loss evolutionary tracks 
were originally published in a series of papers by the 
Geneva group (Schaller et al. 1992; Schaerer et al. 1993a, 
1993b; Charbonnel et al. 1993). These models adopt 
mass loss rates throughout the HR diagram from de Jager 
et al. (1988) that are scaled with metallicity according 
to the models of Kudritzki et al. (1989), where M oc 
^0-5. Wolf-Rayet (WR) stars are an exception - the mass 
loss rates for these stars are taken from Langer (1989) 
and Conti (1988) and include no correction for initial 
metallicity effects. 
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The "high" mass loss evolutionary tracks (HIGH) , pub- 
lished in Mcynct et al. (1994), include enhanced mass 
loss rates, meant to more accurately reproduce observa- 
tions of low-luminosity Wolf-Rayct stars and blue-to-red 
supergiant ratios in the Magellanic Clouds (Schaller et 
al. 1992, Meynet 1993). The adopted mass loss rates 
are derived by doubling the rates adopted by the "stan- 
dard" grid from dc Jagcr et al. (1988), as well as dou- 
bling the "standard" mass loss rate assumed for late-type 
WN-type WR stars. The mass-loss rates for early-type 
WN WR stars and later stages of WR stars (WC and 
WO) were left unchanged - for a complete discussion, see 
Meynet et al. (1994). Again, the mass loss rates of WR 
stars are left unc:orrcc:tccl for initial mctallicity effects. 

While many advances have since been made in our un- 
derstanding of stellar physics, adopting these tracks in 
our stellar population synthesis models is still scientif- 
ically sound. The STD mass loss tracks are the more 
applicable of the two when considering recent work on 
the effect that wind clumping has on mass loss rates 
(Crowther et al. 2002); however, the HIGH mass loss 
tracks produce a reasonable approximation of the en- 
hanced mass loss rates resulting from the effects of ro- 
tation, when surface mixing results in an earlier start of 
the WR phase (Meynet, private communication). Rota- 
tion is an important component of stellar evolution that 
is expected to have considerable influence on the Star- 
burstQQ ionizing spectrum and the agreement of these 
models with observations at low metallicity in particular 
(Leitherer 2008). Since the effects of rotation are not 
explicitly included in these tracks, the HIGH mass loss 
tracks can be considered a more appropriate approxima- 
tion of the rotation-driven mass loss undergone by mas- 
sive stars. 

2.3. Starburst99 Ionizing Spectra 

The far-ultraviolet (FUV) ionizing spectra produced 
by Starburst99 are primarily influenced by age and 
metallicity. The effects of a changing model age derive 
largely from the evolution of the massive stellar popula- 
tion, and are easily examined in our models which adopt 
an instantaneous burst model of star formation - this 
star formation history allows us to observe the effects of 
a single stellar population that is formed at Myr and 
evolves uniformly. 

By contrast, for the continuous SFH models, we set 
the age constant at 5 Myr. This age describes an active 
(emission-line) star-forming galaxy, where the number of 
stars being formed is equal to the number of stars dying. 
At younger ages this equilibrium is not yet reached, and 
at older ages there is little to no evolution in the FUV 
ionizing spectrum produced by the stellar population. 
This is consistent with the evolution of continuous SFH 
models for starburst galaxies, as described in Kewley et 
al. (2001). 

Figure 1 (left) shows the FUV spectra generated when 
adopting an instantaneous burst SFH and the HIGH evo- 
lutionary tracks, plotted from Myr to 6 Myr in 1 Myr 
increments for our full range of metallicities. The hard- 
ness of the spectra decreases with age for this star forma- 
tion history, most noticeably in the higher-energy regime 
of the spectra (100-300 A), a result of the massive stellar 
population evolving out of the OB main sequence phase. 

It is apparent that the behavior of the FUV spec- 



trum differs dramatically across the different metallici- 
ties. The low-metallicity spectra maintain a significantly 
harder ionizing spectrum throughout their evolution as 
compared to the Z = Zq and Z = 2Zq model spec- 
tra. Since high-metallicity stars spend a larger fraction 
of their high-energy photons ionizing their own atmo- 
spheric metals (chic to the increased effects of line blan- 
keting), there is a resulting depletion of high-energy pho- 
tons available to ionize the surrounding ISM, which leads 
to a softer radiation field (Snijders et al. 2007). The effec- 
tive temperatures of massive stars are also higher at lower 
metallicities across similar spectral types due to a shift of 
the evolutionary limits of hydrodynamic equilibrium (the 
Hayashi limit) to warmer temperatures at lower metal- 
licities (Elias et al. 1985; Levesque et al. 2006), resulting 
in a hotter environment and harder spectrum in the case 
of our low-metallicity spectra. Finally, main sequence 
lifetimes are longer at low metallicities as a result of 
lower mass loss rates, leading to a greater amount of time 
spent in the hot main-sequence evolutionary phases and 
a larger contribution to the ionizing radiation field that 
extends to later ages (e.g. Meynet et al. 1994, Maeder & 
Conti 1994). 

Figure 1 (left) also indicates that the Wolf-Rayet (W- 
R) phase for these galaxies contributes to a hardening of 
the ionizing spectrum, producing a distinctive bump in 
the high-energy region of the spectrum. The net effect of 
this population is only visible in the high-Z ionizing spec- 
tra, appearing from 3 to 5 Myr in the Z = Zq spectrum, 
and from 4 to 6 Myr in the Z = 2Zq. The wind-driven 
ISM enhancement by W-R stars is stronger and longer 
at high metallicities, and the minimum mass required 
to reach the W-R stage decreases at higher metallicities 
(taken to be 85 Mq a,t Z = 0.05^0 as compared to 4OM0 
at Z = Zq for the Geneva models used here; Schaller et 
al. 1992). 

However, the most pronounced change in the spectra 
occurs at ages later than 5 Myr. For the 6 Myr model, 
there is a dramatic drop in the hardness of the spectra 
at ~ 225 A, coinciding almost perfectly with the ion- 
izing wavelength of [OIII]. This behavior is consistent 
with the origin of the FUV spectrum for this SFH. An 
instantaneous burst of star-formation at Myr will re- 
sult in a single coeval population of massive stars as the 
sole source of ionizing radiation, decreasing the hard- 
ness of the FUV spectrum at later ages as the hottest 
massive stars (and therefore the dominant contributors 
of ionizing radiation), evolve rapidly off the main se- 
quence and end their lives. The lifetimes employed in 
the Geneva models give the explanation for the behav- 
ior of the FUV spectra shown here; the main-sequence 
(H-burning phase) lifetimes of a typical hot massive star 
of initial mass 40 Mq all terminate at an age of 5 Myr 
or earlier (Schaller et al. 1992, Schaerer et al. 1993b, 
Charbonnel et al. 1993). As a result, beyond 5 Myr 
the dominant massive stellar population is comprised 
mainly of lower-mass (< 25 Mq) stars that are beginning 
to evolve redwards on the Hertzsprung-Russell diagram, 
and producing insufficient ionizing radiation at shorter 
wav clcngths (< 225 A) to generate the line fluxes we ob- 
serve for highly ionized species. Because of this lack of 
ionizing photons, we restrict our photoionization grids 
generated for an instantaneous star formation history to 
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5 Myr and younger. 

In practice, stellar populations in star-forming galax- 
ies are thought to originate from episodic bursts of star 
formation, or from a continuous SFH. For example, Izo- 
tov & Thuan (2004) find evidence of current ongoing star 
formation with an age of about 4 Myr in the extremely 
low-mctallicity galaxy I Zw 18. In addition, they de- 
tect supergiant populations that indicate an intense star 
formation episode occuring 10-15 Myr ago and evidence 
of still older stellar populations with ages of hundreds 
of Myr. Noeske et al. (2007a, 2007b) examine 2905 
star-forming galaxies from the All Wavelength Extended 
Groth Strip International Survey (AEGIS), and find that 
the dominant star formation mechanism in these galaxies 
appears to be a gradual decline of the average star for- 
mation rate, as opposed to a series of episodic bursts that 
decrease in frequency. Previous grids of stellar popula- 
tion synthesis models also find that the assumption of a 
continuous SFH, with the presence of underlying older 
stellar populations, is more appropriate than a single 
instantaneous burst in most cases (Kewley et al. 2001, 
Moy et al. 2001, Fernandes et al. 2003, Barton et al. 
2003). This suggests that, rather than modeling a single 
burst of star- formation at Myr (the current approach in 
Starburst99 for a burst-like SFH), a continuous SFH or 
multiple bursts of star formation through time would be 
a more realistic treatment for modeling spectra of star- 
forming galaxies (see also Lee et al. 2009). The instanta- 
neous burst models presented here can therefore be em- 
ployed mainly as an approximation of a burst-like SFH, 
representing a lower limit on star formation, while our 
continuous star formation models can be considered an 
upper limit. 

After examining the behavior of the FUV spectra when 
adopting the HIGH evolutionary tracks, we can then 
compare these to the FUV spectra produced from the 
STD evolutionary tracks, subtracting the two spectra 
to examine their differences in greater detail (Figure 1, 
right). In this comparison, we can see that the largest 
differences occur at higher metallicities {Z = OAZq, Zq, 
and 2Zq) and in the 3 to 6 Myr age range, coincid- 
ing with the ionizing flux contribution of the mass-loss- 
sensitive Wolf-Rayet phase. Additional discrepancies at 
later ages focus on the age of the 225A drop, which is 
also expected to vary slightly with mass loss. 

Finally, we find that the FUV spectra shown here are 
harder than those shown in Kewley et al. (2001) for both 
the PEGASE and the Starburst99 stellar population syn- 
thesis models. The increased hardness between 225A 
and lOOOA is particularly notable, with log(AFA, ergs 
s""'^ Mq'^) « 25 to 30 for the solar-metallicity Kewley et 

al. (2001) models as compared to log(AFA, ergs s~^ Mq^) 
« 35 to 45 for the solar-metallicity models shown here. 
This corresponds to the 1 to 4 Ryd region that Kewley 
et al. (2001) cite as potentially benefiting from more rig- 
orous treatments of continuum metal opacities, which we 
have adopted here via the inclusion of the Pauldrach et 
al. (2001) and Hillier & Miller (1998) model atmospheres. 
These FUV spectra are therefore a notable improvement 
over the Kewley et al. (2001) grids, with higher ionizing 
fluxes expected to more accurately reproduce the emis- 
sion line fluxes observed in star-forming galaxies. 
Figure 2 shows the relative ionization fractions for a 



number of species produced by Mappings HI, plotted as 
a function of relative distance from inner surface of the 
nebula. This figure illustrates that the [SII] fiux in par- 
ticular is a good tracer of the hardness of the ionizing 
radiation field. This is largely due to the effects of [SHI] 
ionization. In a soft ionizing radiation field [SHI] will be 
ionized out to a greater distance in the nebula, decreas- 
ing the ionization fraction of [SII] as a result. Conversely, 
in a hard ionizing radiation field, [SHI] will be ionized 
very close to the inner surface of the nebula and [SII] 
will dominate the ionization fraction. As a result, [SIl] 
can be used as a powerful diagnostic tool in these mod- 
els, particularly as it is a commonly-detected emission 
feature in star-forming galaxy spectra. 

3. OPTICAL EMISSION LINE DIAGNOSTICS 

The final output of the Starburst99/Mappings HI code 
is a model galaxy emission spectrum. The model spectra 
calculated for our grid of input parameters allows us 
to probe the ISM properties of these galaxies through 
the use of emission line diagnostics. In this work we 
employ the [NII]A6584/Ha, [NII]A6584/[OII]A3727, 
OIII]A5007/H/3, [OIII]A5007/[OII]A3727, and 
SII]AA6717,6730/Ha line ratios. As these fluxes 
can vary with the age of the young stellar population, 
we examine the evolution of each ratio as a function of 
time across our instantaneous burst SFH model grids, 
and consider this evolution in the context of the FUV 
ionizing spectra shown in Figure 1. For these ratios we 
assume an ionization parameter g = 2 x 10* cm s~^; 
this corresponds to log U ~ —2.2, our closest ionization 
parameter to the log U = —2.3 value adopted by Rigby 
& Rieke (2004) in their models. 

3.1. [Nil] /Ha 

The [Nil] A6594/Ha ratio correlates strongly with both 
metallicity and ionization parameter, and is useful in di- 
agnostic diagrams comparing these parameters (Veilleux 
& Osterbrock 1987, Kewley et al. 2001). The [Nil] 
flux at low metallicities is dominated by the abundance 
of primary nitrogen (Chiappini et al. 2005, Mallery et 
al. 2007). Primary nitrogen production is largely in- 
dependent of metallicity and occurs predominantly in 
intermediate-mass stars (Matteucci & Tosi 1985, Mat- 
teucci 1986); however, some primary nitrogen produc- 
tion by massive stars is thought to be dependent on stel- 
lar mass and metallicity (Chiappini et al. 2005, 2006; 
Mallery et al. 2007). At higher metallicities the produc- 
tion of secondary nitrogen becomes prevalent (AUoin et 
al. 1979, Considere et al. 2000, Mallery et al. 2007); sec- 
ondary nitrogen is synthesized from carbon and oxygen 
originally present in the star and is therefore proportional 
to abundance. At low metallicities, nitrogen abundance 
is also relatively sensitive to the star formation history of 
the galaxy, introducing some scatter into the relation be- 
tween [Nil] /Ha and metallicity (Kewley & Dopita 2002). 
The relatively low ionization potential of [Nil], allowing 
production of this line at the lower-energy end of the 
FUV ionizing spectrum, also makes this ratio sensitive 
to ionization parameter (Kewley & Dopita 2002). 

The evolution of this ratio with age is shown in Figure 
3a, where we plot all five metallicities and the HIGH and 
STD mass loss rates for the instantaneous burst SFH 
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models. We can see an increase of the [Nil] /Ha ra- 
tio with metaUicity. At the lowest mctallicities {Z — 
0.05Zq,Z = 0.2Zq) the [NIIJ/Ha ratio increases by 
~80%-160% at later ages (> 3 Myr), while at higher 
metallicities the ratio varies across all ages, showing de- 
viations from the mean of up to ~75%-85%. We can also 
see that [Nil] /Ha becomes double- valued at later ages 
(> 4 Myr) for the highest metallicities {Z — Zq and 
Z = 2Zq). To understand the variation of this ratio with 
age for an instantaneous burst SFH, we can consider the 
FUV ionizing spectrum. In the instantaneous burst FUV 
spectra, we sec a larger variation with age at the ionizing 
threshold of [Nil] as compared to the continuous burst 
FUV spectra, a conscqiicncc of the varying contribution 
from the coeval massive star population. This variation 
also becomes more prominent at higher metallicities, ex- 
plaining the increased variation of the [Nil] /Ha ratio in 
that regime. 

3.2. [NII]/[OII] 

The [Nil] A6594 and [OH] A3727 hncs have very similar 
ionization potentials, making this ratio almost indepen- 
dent of the ionizing parameter of the radiation field. This 
ratio also avoids the pitfall of being double- valued with 
abundance, scaling smoothly from high to low abundance 
(van Zee et al. 1998) and rendering it a reliable means 
of isolating metaUicity on diagnostic diagrams (Dopita et 
al. 2000). [Nil] /[Oil] correlates very strongly with metal- 
licity above Z = OAZq; [NII]'s status as a secondary 
element causes it to more strongly scale with increas- 
ing metaUicity than [Oil] . In addition, the lower electron 
temperature at high metaUicity produces fewer thermal 
electrons, decreasing the number of coUisional excitations 
in the high-cnergy-thrcshold [OH] feature as compared to 
the lower-energy [Nil] features (Kewley & Dopita 2002). 
At lower metallicities {Z < 0.23^©) this ratio's utility as 
a tracer of metaUicity decreases slightly, since nitrogen 
and oxygen are both primary nucleosynthesis elements 
and begin to scale uniformly with metaUicity as a result 
(Dopita et al. 2000). 

The evolution of this ratio is shown in Figure 3b. This 
ratio shows increased variation with age at higher metal- 
licities, similar to the behavior observed in the [Nil] /Ha 
ratio. As described in the [Nil] /Ha case, there is a vari- 
ation with age in the FUV spectrum at the [Nil] ionizing 
threshold that increases at higher metallicities, leading 
to increased variations in the flux of the [Nil] line. In 
this case, we see the same behavior in the FUV spec- 
trum at the ionizing threshold of [Oil] . The slope of the 
FUV ionizing spectrum between these two thresholds in- 
creases with metaUicity and age. This change in the rela- 
tive ionizing flux between the two species leads to greater 
variation in the [Nil] /[Oil] flux at higher ages and higher 
metallicities. We can, however, see that the [Nil] /[OH] 
ratio increases significantly with metaUicity, demonstrat- 
ing that on a diagnostic diagram a clear separation with 
metaUicity will be evident across the [Nil] /[Oil] axis. 

3.3. [0IIIJ/HI3 

The [OIII] A5007/H/3 ratio is sensitive to the hardness 
of the ionizing radiation field, and a useful means of trac- 
ing the ionizing parameter of a galaxy (Baldwin et al. 
1981). However, this ratio is sensitive to metaUicity and 



double-valued with respect to abundance, although it is 
far more sensitive to ionization parameter at sub-solar 
metallicities (Kewley et al. 2004). The close proximity 
of these emission lines also renders this ratio relatively 
insensitive to reddening corrections. 

The evolution of the [0III]/H/3 ratio is shown in Fig- 
ure 3c. We can see immediately that, as described in 
Kewley et al. (2004), [0III]/H/3 is double- valued with 
metaUicity, reaching its highest value at Z = OAZq and 
decreasing for both lower and higher metallicities. We 
can also see that the value of the ratio decreases with 
age. This decrease is gradual at the lowest metallicities 
but is interrupted by local maxima at 3 to 5 Myr as 
metaUicity increases. This age range conibiucid with the 
behavior of the FUV ionizing spectrum indicates that 
these local maxima in the [0III]/H/3 ratio are due to the 
short-lived contribution of the Wolf-Rayet phase to the 
flux of the [OIII] line. We can see sX Z = Q.2Zq that 
this contribution from the Wolf-Rayet popiUation is ob- 
served in the HIGH mass loss models but not the STD 
models; recall that at low metallicities in particular the 
HIGH mass loss models are expected to produce an ear- 
lier Wolf-Rayet population than the STD models. The 
overall decrease with age in the [0III]/H/3 flux can be at- 
tributed to the gradual decrease in hardness of the FUV 
ionizing spectrum at the ionization threshold of [OIII]. 

3.4. [OIII]/[OII] 

The [OIII] A5007/[OII] A3727 serves a similar function 
as the [0III]/H/3 ratio, and is a commonly-employed ion- 
ization parameter diagnostic. Dopita et al. (2000) and 

Kewley & Dopita (2002) do find that the [OIII] /[Oil] is 
sensitive to metaUicity as well, with higher abundances 
corresponding to lower values, but both the metaUicity 
and ionization parameter relations are monotonic and do 
not result in this ratio being double- valued. 

The evolution of the [OIII] /[OH] ratio is shown in Fig- 
ure 3d, and supports the conclusions of Dopita et al. 
(2000); we can sec that the value of this diagnostic ra- 
tio decreases as metaUicity increases. The evolution of 
this diagnostic with age is quite similar to that of the 
[0III]/H/3 ratio. As in the case of [OIII]/H/?, the [OIII] 
flux here is affected by the slowly decreasing hardness of 
the FUV spectra and the contributions of the Wolf-Rayet 
phase. By contrast, the flux of the FUV spectrum stays 
relatively constant at the ionizing threshold of [Oil] for 
both SFHs, isolating the behavior of the [OIII] flux in 
the evolution of this ratio. 

3.5. [SIIJ/Ha 

The [SII] (A6716 -|- A6731)/Ha line ratio is a useful 

means of tracing the hardness of the photoionizing spec- 
trum present in a galaxy. While the ratio retains a slight 
dependence on the ionization parameter, we define this 
value at the innermost boundary of the nebula. [SII] , by 
contrast, forms at a greater distance from this boundary 
in a partially ionized zone. As the hardness of the ion- 
izing radiation field increases the partially ionized zone 
becomes more extended, generating increased [SII] flux. 

The evolution of the [SII] /Ha ratio is shown in Figure 
3e. We can see that the value of this ratio increases 
gradually with abundance across the lower metal- 
licities but becomes degenerate at high metallicities 
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{Z ^ Zq,Z = 2Zq). We see that the [SII]/Ha flux 
increases with age and undergoes large fluctuations at 
higher metallicities, showing a ~90%-95% deviation 
from the mean. This is, once again, a signature of the 
increased variation with metallicity seen in the FUV 
ionizing spectrum at the [SIT] ionizing threshold. Out 
of all of oiu- diagnostic ratios, the [SII] flux is also the 
most sensitive to small variations in the hardness of the 
FUV ionizing spectrum. 

We find that all the ratios are degenerate with 
age, particularly at high metallicity. It is therefore 
important to determine the young stellar population 
age for a galaxy before drawing any robust conclusions 
through the use of the instantaneous burst SFH models. 
We also find a measurable difference between the HIGH 
and STD mass loss rate models. This once again 
demonstrates that an instantaneous burst SFH traces 
the effect that a single coeval population of massive stars 
has on the FUV ionizing spectrum and clearly illustrates 
differences in the emission line spectrum imposed by the 
different stellar evolutionary tracks. 

4. EMISSION LINE DIAGNOSTIC DIAGRAMS 

We wish to test our Starburst99/Mappings models 
against observed spectra from a variety of galaxy pop- 
ulations. The local {z < 0.1) galaxies we selected for our 
comparison samples are described below. 

Sloan Digital Sky Survey (SDSS) galaxies: To com- 
pare our models to the general local galaxy population, 
wc plot emission line ratios determined for a sample of 
emission-line galaxies from SDSS described in Kewley 
et al. (2006). These galaxies were originally taken from 
Data Release 4 of SDSS (Adelman-McCarthy et al. 2006) 
and restricted to 85224 galaxies with a S/N > 8 in the 
strong emission lines (Tremonti et al. 2004) and a redshift 
range between 0.04 < z < 0.1. The lower limit of this 
redshift range corresponds to an aperture covering frac- 
tion of 20%, the minimum required to avoid domination 
of the spectrum by aperture effects (Kewley et al. 2005, 
Kewley & Ellison 2008) . The spectra were acquired with 
3" diameter fibers. We further restrict this sample to 
the 60920 galaxies classified as star-forming by Kewley 
et al. (2006). This classification is based on the criteria 
of Kauffmann et al. (2003) derived from the [Nil] /Ha 
vs. [0III]/H/3 diagnostic diagram, and the Kewley et al. 
(2001) criteria based on the [SII]/Ha vs. [0III]/H/3 and 
[OI]/Ha vs. [OIIIJ/H/3 diagnostic diagrams; equations 
are given in Kewley et al. (2006). This classification re- 
moves contaminations from composite galaxies, Seyferts, 
LINERs, and ambiguous galaxies that cannot be defini- 
tively classified based on the Kauffmann et al. (2003) and 
Kewley et al. (2001) criteria. 

Nearby Field Galaxy Survey (NFGS) galaxies: To en- 
sure that the local galaxy population comparison is not 
affected by any residual aperture effects present in the 
SDSS sample, we include a sample of galaxies from 
NFGS. NFGS has integrated emission-line spectropho- 
tometry of 196 galaxies that span from —14 < Mb < —22 
and include the full range of Hubble sequence morpholo- 
gies, with a median redshift of z — 0.01 and a maxi- 
mum redshift of ^ ~ 0.07 (Jansen et al. 2000a, 2000b). 
The galaxies were selected from the CfA redshift cata- 
log, which has a limiting blue photographic magnitude 



of mz — 14.5 (Huchra et al. 1983) and observed with 
the FAST spectrograph at the Whipple 1.5 m telescope 
(Jansen et al. 2000b). We limit this sample to 95 star- 
forming galaxies by applying the emission-line criteria for 
star-forming galaxies in Kewley et al. (2006) described 
above for the SDSS sample. 

Blue Compact Galaxies (BCGs): To compare our mod- 
els to a sample of galaxies that show evidence of burst- 
like star formation histories, we include a sample of blue 
compact galaxies (BCGs) from Kong & Cheng (2002) 
and Kong et al. (2002) . Their sample consists of 97 galax- 
ies selected from the of Gordon & Gottesman (1981), 
Thuan & Martin (1981), Kinney et al. (1993), and Thuan 
et al. (1999) and limited to Mb < —17 mag. The spectra 
were observed with the OMR spectrograph on the 2.16m 
telescope of the Beijing Astronomical Observatory. Here 
we restrict our sample to those galaxies with a complete 
set of spectral line fiuxes that satisfy the star-forming 
galaxy emission line criteria in Kewley et al. (2006), giv- 
ing us 36 star- forming BCGs in our final sample. Note 
that, for consistency with our other star-forming galaxy 
samples, we use the Kewley et al. (2006) criteria rather 
than the classifications used in Kong et al. (2002) to 
remove AGN and non-emission-line galaxies from their 
sample. This sample of BCGs has a redshift range of 
0.003 < z < 0.029, with a median redshift of z = 0.016. 

Metal-Poor Galaxies (MPGs): Finally, to examine the 
agreement of our stellar population synthesis and pho- 
toionization models with low metallicity galaxies, we in- 
clude a sample of 10 MPG spectra from Brown et al. 
(2008). The 10 MPG spectra are part of a survey of 38 
MPG candidates, selected from Data Release 4 of SDSS 
based on restrictions on their [u' — g')o, [g' — r')o, and 
{r' — i')o colors, a magnitude limit of g' < 20.5, and re- 
moval of HII regions in nearby galaxies through visual in- 
spection. The remaining MPG candidates were observed 
with the Blue Channel spectrograph at the MMT tele- 
scope. Brown et al. (2008) find 10 MPGs in this sample, 
with log(0/H) -h 12 < 8 based on their electron temper- 
ature metallicities (Izotov et al. 2006). These 10 galaxies 
have a redshift range of 0.03 < z < 0.081 with a median 
redshift of z = 0.073. 

All of the emission line fluxes used in our comparison 
samples have been corrected for total intrinsic attenu- 
ation by dust in the direction of the galaxy. This was 
based on the Ha/H/3 emission line ratio, assuming the 
Balmer decrement for case B recombination (Ha/H/? = 
2.87 for T = 10'' K and - 10^ - 10"* cm-'\ following 
Osterbrock 1989) and the Cardelh et al. (1989) reddening 
law with the standard total-to-selective extinction ratio 
Rv = 3.1. 

These comparisons samples allow us to examine the 
agreement of our stellar population synthesis and pho- 
toionizaion models when applied to the local galaxy pop- 
ulation, including nearby low-metallicity galaxies. We 
plot these comparison samples along with our models on 
a series of optical emission line ratio diagnostic diagrams, 
utilizing the ratios detailed in Section 3. 

4.1. [Nil] /Ha vs. [GUI] /HP 

In Section 3, we note that [Nil] /Ha correlates strongly 
with metallicity as well as with ionization parameter 
(VeiUeux & Osterbrock 1987, Kewley et al. 2001; see 
Section 3.1), while [OIII]/H/J is primarily a measure of 
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ionization parameter with a degenerate dependence on 
metallicity (Baldwin et al. 1981, Kewlcy et al. 2004; sec 
Section 3.3). The [NII]/Ha vs. [0III]/H/3 grid of Bald- 
win et al. (1981) allows us to examine these two ISM 
properties by examining the evolution of the grids with 
age and comparing the agreement of our models to ob- 
served galaxy spectra. In Figure 4 we plot the evolution 
of the [Nil] /Ha vs. [OIIIJ/II/3 diagnostic grid for our 
instantaneous burst SFH models ranging from 0.0 to 5.0 
Myr, while Figure 5 shows the diagnostic diagram assum- 
ing a continuous SFH and an age of 5.0 Myr; both figures 
compare the model grid to the SDSS, NFGS, BCG, and 
MPG samples described above. 

Following the observed evolution of these diagnostic 
ratios in Section 3, we can see that the instantaneous 
burst models in Figure 4 show a substantial decrease in 
the [0III]/H/3 flux with age, attributable to the signifl- 
cant decrease in the hardness of the FUV spectrum with 
age. When c;omparing these grids to our galaxy sample, 
we see satisfactory agreement for the youngest instan- 
taneous burst models (0.0 Myr, 1.0 Myr); these mod- 
els accommodate 69% of the SDSS galaxies, 78% of the 
NFGS sample, 58% of the BCGs, and 86% of the MPGs, 
a great improvement over previous model grids that show 
poor agreement with metal-poor galaxies. However, this 
agreement rapidly deteriorates with age as the [OIIIJ/H/3 
ratio in the models decreases and the fluctuating contri- 
butions of the Wolf-Rayet phase are introduced at 3.0 
to 5.0 Myr. The MPG sample in particular is in poor 
agreement with the models at 3.0 to 5.0 Myr. At later 
ages (4.0 Myr, 5.0 Myr) the HIGH mass loss models show 
slightly better agreement with the data as compared to 
the STD models. By comparison, we find that the contin- 
uous SFH models in Figure 5 show a better agreement 
with the data at 5.0 Myr, although deficiencies in the 
model [OIIIJ/H/9 fluxes are apparent. 

From this diagnostic, it appears that our sample popu- 
lation is restricted to the higher metallicities included in 
our models (with the MPGs and a small sample of out- 
liers as the only exception) and spans the full range of 
ionization parameters included in our models, again with 
only a small sample of outliers. However, the double- 
valued nature of this diagnostic makes it an impractical 
sole means of drawing strong conclusions about galaxies' 
metallicities and ionization parameters. 

4.2. [NII]/[OII] vs. [OIII]/[OII] 

In Figure 6 we plot the evolution of [Nil] /[OH] vs. 
[QUI] /[on] with age for the instantaneous burst SFH 
models and compare them to our galaxy samples; in Fig- 
ure 7 we again plot the 5.0 Myr continuous SFH models. 
We can see from the model grids plotted in Figure 6 and 
Figure 7 that [NII/OII] is primarily sensitive to metallic- 
ity, while [OIII/OII] is primarily sensitive to ionization 
parameter, following from discussion of these diagnostic 
ratio in Sections 3.2 and 3.4. 

This diagnostic shows much more consistent agreement 
with our samples - the best-fit model grids at O.OMyr ac- 
commodate 91% of the SDSS galaxies, 66% of the NFGS 
sample, 86% of the BCGs, and 86% of the MPGs. For 
the instantaneous burst models, we see a rapid decrease 
with age in the [QUI] /[OH] ratio and fluctuations from 
the Wolf-Rayet phase contribute at 3.0 to 5.0 Myr. The 
1.0 to 2.0 Myr instantaneous burst models also appear to 



be somewhat deficient at higher metallicities in the case 
of the SDSS sample, while the 5.0 Myr grid shows poor 
agreement with the low-metallicity MPG sample. At 5.0 
Myr, the HIGH mass loss models are also found to be 
in better agreement with our galaxy samples than the 
STD models. By contrast, the continuous SFH models 
show consistent agreement (66%)-95%i) with o\iv galaxy 
samples at 5.0 Myr for both the HIGH and STD mass 
loss evolutionary track models, including all but the most 
metal-poor galaxy in the MPG sample. 

In this non-degenerate diagnostic we can see that the 
SDSS galaxies span the full range of metallicities adopted 
in our models, although they appear to be restricted to 
the lower ionization parameters ((? < 8 x 10^ cm s^^). 
The NFGS and BCG samples span a similar range in ion- 
ization parameter, but do not extend to the twice-solar 
metallicity regime of our models. Interestingly, the MPG 
sample appear to coincide with both the lowest metallic- 
ities [Z = 0.05^0 to Z = 0.2^0) and the highest ioniza- 
tion parameters present across our comparison samples, 
along with several of the most metal-poor BCGs, ranging 
from q = 8 X lO*" cm s^^ to q — 2 x 10^ cm s~-^. 

4.3. [SII]/Ha vs. [0III]/Hf3 

In Figure 8 we plot the evolution of the [SII]/Hq! vs. 
[0III]/H/3 diagnostic ratio for our models ranging from 
0.0 to 5.0 Myr assuming an instantaneous burst SFH, 
while Figure 9 shows the diagnostic at 5.0 Myr for a 
continuous SFH. In both figures we compare the model 
grid to our galaxy samples. 

Both the instantaneous burst and continuous SFH 
models show partial agreement with the galaxy samples; 
however, it is clear in these diagrams that the grid is not 
able to accommodate a large number of these galaxies. 
For the best-fit instantaneous burst models at 0.0 Myr, 
the models only show agreement with 61% of the SDSS 
galaxies, 26% of the NFGS galaxies, 33% of the BCGs, 
and 50% of the MPGs. Agreement is better across the 
[SII]/Ha axis, but due to the double- valued nature of the 
diagram we cannot isolate the behavior of this line ratio. 
In the instantaneous burst models, we see fluctuations 
in the [0III]/H/3 ratio with age and progressively poorer 
agreement with the galaxy samples, beginning with a 
failure to accommodate any of the galaxies in the MPG 
sample by 3.0 Myr but progressing to only agreeing with 
35% of the SDSS galaxies, 13% of the BCGs, and 6% 
of the NFGS galaxies by 5.0 Myr for the HIGH mass 
loss grid; the STD mass loss grid shows a ~0% agree- 
ment with our galaxy samples at 5.0 Myr. The continu- 
ous SFH models maintain a somewhat more satisfactory 
agreement at 5.0 Myr and show little distinction between 
the HIGH and STD mass loss rates, but the grids are still 
insufficient and show on 30% agreement with the MPG 
sample. As discussed in Section 3.5, use of the [SII]/Ha 
diagnostic allows us to examine the hardness of the FUV 
ionizing spectrum; Figures 8 and 9 suggest that a gen- 
erally harder FUV ionizing spectrum is required from 
the models to produce a more extended partially ionized 
zone in the theoretical nebula. 

4.4. Comparison with Previous Model Grids 

The stellar population synthesis and photoionization 
models presented in Dopita et al. (2006) are compara- 
ble in many ways to the models presented in this work. 
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They use the Starburst99 stellar population synthesis 
code, adopting the Pauldrach et al. (2001) and Hillier 
& Miller (1998) model atmospheres that include treat- 
ments of metal opacities and the five metallicities of the 
Geneva HIGH evolutionary tracks. They also use the 
latest generation of Mappings III to compute their pho- 
toionization models. The Dopita ct al. (2006) models 
adopt a spherical geometry in the photoionization mod- 
els, while this work assumes a plane-parallel nebular ge- 
ometry - these different Mappings geometries have been 
compared and produce equivalent results (Kewley et al. 
2001). There are, however, two noteworthy differences 
between the models of Dopita et al. (2006) and those 
presented here. 

First, the models of Dopita et al. (2006) do not take 
ionization parameter (g) as one of their free parameters; 
instead, q is replace by TZ, a parameter representing the 
ratio of the mass of the central aging star cluster to the 
pressure of the surrounding ISM. This allows q to vary 
with age (as well as with metallicity) . In this treatment, 
TZ is fixed at a variety of values (—6, —4, —2, 0, and 2) and 
a sequence of model HII regions are computed at each of 
the Geneva metallicities with ages that are increased in 
increments of 0.5Myr, up to a maximum age of 6.5Myr. 

Second, Dopita et al. (2006) generate models of indi- 
vidual HII regions. To model the spectra of star- forming 
galaxies, they integrate the fluxes of a model HII region 
for each of their model ages, essentially considering star- 
forming galaxies to have spectra that consist of contri- 
butions from multiple HII regions at different ages. This 
removes age as a free parameter in these models, restrict- 
ing the free parameters to metallicity and TZ (similar to 
our models with free parameters of metallicity and q). 

As compared to the Kewley et al. (2001) models, the 
primary difference in this work lies in the choice of model 
atmospheres. While this work adopts the the Pauldrach 
et al. (2001) and Hillier & Miller (1998) model atmo- 
spheres, the Kewley et al. (2001) grids utilize the Leje- 
une et al. (1997) and Schmutz et al. (1992) models, which 
do not include treatments of metal opacities. Kewley et 
al. (2001) cite this as a shortcoming of their models and 
suggest that the inclusion of continuum metal opacities 
will result in a harder FUV spectrum, a prediction that 
is supported by our models (see Section 3). The Kewley 
et al. (2001) models also include slightly different treat- 
ments of rig (rie = 350 cm~^ as compared to the rig = 100 
cm^^ models shown here) and age (adopting a continu- 
ous SFH model age of 8.0 Myr as opposed to the 5.0 Myr 
used in this work). However, both of these differences are 
expected to have negligible effects on the diagnostic grids 
produced by these models. 

In Figure 10 (left) we compare our 0.0 Myr instanta- 
neous burst SFH models to the Kewley et al. (2001) and 
Dopita et al. (2006) models on the emission line diagnos- 
tic diagrams described above, considering the agreement 
of all three grids with the SDSS and MPG galaxy sam- 
ples. The Kewley et al. (2001) models also assume a 0.0 
Myr instantaneous burst SFH and range from q = Ix 10^ 
to g 3 X 10^; the Dopita et al. (2006) models assume an 
instantaneous burst SFH and range from 0.0 Myr to 4.0 
Myr {q is not a free parameter in the Dopita et al. 2006 
models). For the [Nil] /Ha vs. [0III]/H/3 diagnostic di- 
agram (top left), our models are found to more closely 



agree with the data than either the Kewley et al. (2001) 
or the Dopita et al. (2006) grids. The Dopita et al. (2006) 
do not produce sufficient [OIII] fluxes at higher metallic- 
ities. The Kewley et al. (2001) accommodate all of the 
galaxies but do not properly track the empirical Kauf- 
mann et al. (2003) cut-off for star-forming galaxies that 
has been applied to the SDSS sample, suggesting that 
the [OIII] fluxes produced by these models are unreal- 
istically high. By contrast, the model grids produced 
by this work track the Kaufmann et al. (2003) cut-off 
perfectly. In the case of the [Nil] /[OH] vs. [OIII]/[OII] 
(center left), all three models show agreement with the 
SDSS and MPG samples, although the agreement with 
the MPG sample is better for the Dopita et al. (2006) 
models and this work (6 out of 7) than for the Kewley et 
al. (2001) models (3 out of 7). Finally, the [SII]/Ha vs. 
[0III]/H/3 diagram (bottom left) shows the failure of the 
Dopita et al. (2006) models to produce accurate [SII]/Ha 
ratios for either sample. The Kewley et al. (2001) mod- 
els have higher [OIII] fluxes than this work; however, our 
models extend to higher values of [SII]/Ha in the SDSS 
sample. Across all three of these models, we find that 
changes in the treatment of dust have minimal effects in 
the optical regime; however, this is expected to make a 
more substantial difference in the infrared. 

This comparison is noteworthy when considering that 
the models presented here are found to be an improve- 
ment over, or at least equivalent to, the integrated model 
spectra of Dopita et al. (2006). In our work we are 
modeling galaxy environments in a simple way by as- 
suming the HII regions observed within an aperture can 
be modeled by a luminosity-weighted mean HII region 
represented by a plane parallel model. It is significant 
that this simple approach produces comparable results 
when compared with the more sophisticated treatment of 
modeling the integrated spectra of multiple HII regions 
employed in the Dopita et al. (2006) models. We can 
therefore conclude that treating model galaxies as single 
luminosity- weighted HII regions, particularly when as- 
sume a zero-age instantaneous burst SFH, is a simpler 
and equally effective approach in stellar population syn- 
thesis and photoionization modeling. 

In Figure 10 (right) we compare our 5.0 Myr continu- 
ous SFH models to the Kewley et al. (2001) 8.0 Myr con- 
tinuous SFH models, which again range from q = lx 10^ 
to g 3 X 10^ For the [NHj/Ha vs. [0III]/H/3 diagnos- 
tic diagram (top right), we can see that the agreement 
of our models with the Kaufmann et al. (2003) cut-off in 
the SDSS sample has degraded for these later-age models; 
the Kewley ct al. (2001) maintain higher [OIII] fluxes. In 
the [Nil] /[Oil] vs. [OIII] /[OH] diagnostic diagram (cen- 
ter right), we again see a slightly improved agreement 
with the MPG sample for our models (6 out of 7) as 
compared to the Kewley et al. (2001) models (4 out of 
7); our models appear to extend to lower metallicities for 
the [NII]/[OII] diagnostic ratio. Finally, on the [SII]/Ha 
vs. [0III]/H/3 diagnostic diagram (bottom right) we can 
see that both model grids do a poor job of accommodat- 
ing the SDSS and MPG samples, although the Kewley et 
al. (2001) once again has higher [OIII] fluxes as compared 
to the models in this work. 

This result is surprising considering the speculation by 
Kewley et al. (2001) that the inclusion of model atmo- 
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spheres with metal opaeitics should lead to a harder FUV 
spectrum and strong line fluxes. While we do indeed 
produce a harder FUV spectrum by including the Paul- 
drach et al. (2001) and Hillier & Miller (1998) models in 
our Starburst99 simulations, the effect on the line ratios 
produced by Mappings III appears to be less than an- 
ticipated. Our models do include a more accurate treat- 
ment of physical conditions, and show improvements in 
the agreement with the Kaufmann et al. (2003) crite- 
ria for star-forming galaxies and the metal-poor galax- 
ies along the metallicity-sensitive [Nil] /[Oil] diagnostic 
ratio. However, it is clear that more substantial im- 
provements in models of star-forming galaxies will re- 
quire further systematic changes in existing stellar pop- 
ulation synthesis and photoionization codes. 

5. DISCUSSION AND FUTURE WORK 

We have generated an extensive suite of models for 
star-forming galaxies, utilizing the newest generation of 
the Starburst99 stellar population synthesis code and the 
Mappings III photoionization code. With these codes 
we have constructed a grid of models with a variety of 
metallicities and ionization parameters, adopting both 
an instantaneous burst SFH and a continuous SFH as 
well as both the HIGH and STD mass loss prescriptions 
from the evolutionary tracks of the Geneva group. This 
grids have been made available to the public^, and are 
being used to develop updated metallicity diagnostics as 
part of a current HST theory grant (P.I. Kewley), with a 
publication of new diagnostics expected within the next 
few months. 

We have examined the ionizing spectrum generated by 
Starburst99 for these models, along with the evolution of 
a number of optical emission line diagnostic ratios with 
time. We have selected a number of local (z < 0.1) star- 
forming galaxies as a comparison sample, from Data Re- 
lease 4 of SDSS (Kewley et al. 2006), the NFGS survey 
of Jansen et al. (2000a, 2000b), the blue compact galaxy 
survey of Kong & Cheng (2002) and the MPG galaxy 
sample of Brown et al. (2008). By comparing our mod- 
els to these data on a scries of emission line diagnostic 
diagrams, we are able to make the following conclusions: 

(1) Models that assume a continuous SFH at 5.0 Myr 
produce a harder FUV ionizing spectrum and show bet- 
ter agreement with the observed emission line ratios of 
our galaxy sample as compared to models with an in- 
stantaneous burst SFH at 5.0 Myr. This is in accordance 
with past work that suggests a continuous treatment of 
star formation is more realistic than a single zero-age in- 
stantaneous burst when modeling star-forming galaxies 
(Kewley et al. 2001, Moy et al. 2001, Ferndandes et al. 
2003, Noeskc et al. 2007a, 2007b). 

(2) Assumption of either the HIGH or STD Geneva 
mass loss rates is found to make very little diS'erence in 
the precision of the continuous SFH models; however, 
in the case of the instantaneous burst models, the HIGH 
mass loss rates produce better agreement with our galaxy 
sample, suggesting that an enhanced rate of mass loss is 
more realistic under the assumption of an instantaneous 
burst. 

(3) From the [NII]/Ha vs. [OIII]/H/? and [NII]/[OII] 

^ The models are available for download at 
http : / /www . emlevesque . com/model-grids/ 



VS. [OIII]/[OII] diagnostic diagrams, we find that the 
metallicity and ionization parameter range of our models 
are in agreement with the ISM properties of the galaxies 
in our comparison sample, including our low-metallicity 
galaxy sample. We find that that these metal-poor galax- 
ies appear to include both lower metallicities and higher 
ionization parameters than our general local galaxy sam- 
ple. 

(4) From the [SII]/Ha vs. [0III]/H/3 diagnostic dia- 
gram, it appears that our models still produce an insuf- 
ficiently hard FUV ionizing spectrum that cannot fully 
reproduce some of the line fluxes observed in our galaxy 

sample. 

(5) Our models of a single luminosity-weighted HII 
region are comparable in precision to models such as 
those presented in Dopita et al. (2006), which integrate 
the spectra of multiple HII regions when modeling star- 
forming galaxies. 

It is important to note that these models still include 
several shortcomings that must be considered when ap- 
plying them to star-forming galaxies and considering fu- 
ture work in this area. One ongoing challenge we must 
consider is the required hardness of the FUV ionizing 
spectrum. The Starburst99 stellar population synthesis 
models shown here have a considerably harder FUV ion- 
izing spectrum in the 225A to lOOOA regime than the 
Kewley et al. (2001) models, a regime that is critical 
in ionizing the forbidden optical emission lines used in 
these analyses. This is a consequence of adopting the 
Pauldrach et al. (2001) and HiUier & Miller (1998) model 
atmospheres, which include detailed treatments of metal 
opacities, an improvement originally proposed in Kewley 
et al. (2001). However, our results show that the models 
still do not produce sufHcient flux in the FUV ionizing 
spectrum; this is most evident in the case of the [SII]/Ha 
vs. [OIII]/H/? diagnostic. [SII] in particular requires a 
larger partially ionized zone generated by a harder radia- 
tion field in the models before it can properly reproduce 
the fluxes observed in our galaxy samples. This suggests 
that further systematic changes are required in the stel- 
lar population synthesis models to produce harder FUV 
ionizing spectra. 

One potential means of addressing this issue could be 
the adoption of stellar evolutionary tracks that include 
the effects of rotation on the stellar population, such as 
the new generation of Geneva evolutionary tracks pre- 
sented in Vazquez et al. (2007). Starburst99 outputs gen- 
erated using the z = 0.02 rotating Geneva models have 
been made available to us (Leitherer, personal commu- 
nication), allowing us to examine the effect that rotation 
has on the ionizing radiation fleld produced by this code. 
Figure 11 shows a comparison of the FUV ionizing spec- 
tra generated by Starburst99 when adopting the HIGH, 
STD, and rotating models of the Geneva group; it is ev- 
ident that evolutionary tracks that include rotation pro- 
duce an ionizing spectrum that is in general harder than 
either of the tracks employed in this work, with the differ- 
ence becoming most pronounced at wavelengths shorter 
than ^225A. This is precisely the effect that would be 
anticipated as a result of including rotation in massive 
stellar evolutionary models. As just one example of the 
changes expected with these new tracks, massive stars 
are found to be hotter and more luminous than previ- 
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ously thought, leading to a hardening of the SED that is 
most prominent in the higher-energy regime of the spec- 
trum (Leitherer et al. 2008). In particular, these rota- 
tion effects are expected to become more pronounced at 
lower metallicity, as rotation eventually becomes a domi- 
nant parameter in the evolution of extremely metal-poor 
stars (Hirschi et al. 2008, Leitherer et al. 2008). We plan 
to expand our current grid of evolutionary models by 
generating a new set of Starburst99/Mappings III syn- 
thetic spectra once the full grid of Geneva evolutionary 
tracks with rotation are available. 

In summary, we find that our stellar population syn- 
thesis and photoionization models of a single luminosity- 
weighted HII region, adopting a continuous SFH, repro- 
duce the line ratios of up to 91% of the local population 
of star-formation galaxies, as well as up to 86% of the 
low-metallicity population. These models will allow us 
to study the ISM properties of metal-poor galaxies in un- 
precedented detail, allowing study of the pre-enrichment 
ISM that could shed new light on star formation pro- 
cesses and mechanisms of ISM enrichment (Brown et 
al. 2008). Models of low-metallicity galaxies can also 
be used to probe potential evolutionary mechanisms for 
metal-poor environments (Kewley et al. 2007. Finally, a 



thorough understanding of metal-poor galaxies and their 
stellar populations could prove beneficial to studying the 
host galaxies of long-duration gamma-ray bursts, which 
are thought to be low-metallicity (Stanek et al. 2006, 
Fruchter et al. 2006, Kewley et al. 2007, Modjaz et al. 
2008). These models still include several shortcomings 
that must be considered, in particular their production 
of an insufficiently hard FUV ionizing spectrum that 
specifically affects the [SII] emission line strength. Fu- 
ture model grids that include the effects of rotation on 
the stellar population may help to resolve this issue. 

We would like to thank the anonymous refercx^ for 
extremely helpful and constructive comments regarding 
this manuscript. We gratefully acknowledge our useful 
correspondence with Warren Brown, Margaret Geller, 
Glaus Leitherer, Georges Meynet, Daniel Schaerer, and 
Leonie Snijders. E. Levesque's participation was made 
possible in part by a Ford Foundation Predoctoral Fel- 
lowship. K. Larson's participation in this project was 
made possible in part through the National Science Foun- 
dation's Research Experience for Undergraduates pro- 
gram. L. Kewley and E. Levesque gratefully acknowledge 
support by NSF EARLY CAREER AWARD AST07- 
48559. 



REFERENCES 



Adelman-McCaxthy, J. K. et al. 2006, ApJS, 162, 38 

Allen, M. G., Groves, B. A., Dopita, M. A., Sutherland, R. S., & 

Kewley, L. J. 2008, ApJS, 178, 20 
AUoin, D., CoUin-Souffrin, S., Joly, M., & Vigroux, L. 1979, A&A, 

78, 200 

Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, Pub. A. S. 
P., 93, 5 

Barton Gillespie, E., Geller, M. J., & Kenyon, S. J. 2003, ApJ, 582, 
668 

Bicker, J. & Pritze-v. Alvcnslcbcn, U. 2005, A&A, 443, L19 
Binette, L., Dopita, M. A., & Tuohy, I. R. 1985, ApJ, 297, 476 
Bressan, A., Fagotto, P., Bertelli, G., & Chiosi, C. 1993, A&AS, 

100, 647 

Brown, W., Kewley, L. J., & Geller, M. J. 2008, AJ, 135, 92 
Bruzual, A. G. & Chariot, S. 1993, ApJ, 405, 538 
CardcUi, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Charbonnel, C, Meynet, G., Maeder, A., Schaller, G., Schaerer, D. 

1993, A&AS, 101, 415 
Chariot, S. & Bruzual. G. A. 1991, ApJ, 367, 126 
Chariot, S. & Longhetti, M. 2001, MNRAS, 323, 887 
Chiappini, C, Matteucci, P., Ballerro, S. K. 2005, A&A, 437, 429 
Chiappini, C, Hirschi, R., Meynet, G., Maeder, A., & Matteucci, 

P. 2006, A&A, 449, 27 
Considere, S., Coziol, R., Contini, T., & Davoust, E. 2000, A&A, 

356, 89 

Conti, P. S. 1988. In 0-Stars and WR Stars, NASA SP-497. ed. 

PS Conti, AB Underbill, pp. 81-269. Washington, DC: NASA 
Crowthcr, P. A., Dessart, K., Hillier, D. J., Abbott, D. B., & 

FuUerton, A. W. 2002, A&A, 392, 653 
dc Jager, C, Nieuwenhuijzen, H., van der Hucht, K. A. 1988, 

A&AS, 72, 259 

Dopita, M. A., Kewley, L. J., Heislcr, C. A., & Sutherland, R. S. 

2000, ApJ, 542, 224 
Dopita, M. A. et al. 2006, ApJS, 167, 177 

Elias, J. H., Frogel, J. A., & Humphreys, R. M. 1985, ApJS, 57, 91 
Perland, G. J. 1996, HAZY, a brief introduction to Cloudy, 
University of Kentucky, Department of Physics and Astronomy 

Internal Report 

Perland, G. J., Korista, K. T., Verner, D. A., Ferguson, J. W., 
Kingdon, J. B., Verner, E. M. 1998, PASP, 110, 761 

Fernandes, R. C, Leao, J. R. S., Lacerda, R. R. 2003, MNRAS, 
340, 29 

Fioc, M. & Rocca-Volmerange, B. 1997, A&A, 326, 950 
Fruchter, A. S. et al. 2006, Nature, 441, 463 

Gallagher, J. S., Hunter, D. A., & Bushouse, H. 1989, AJ, 97, 700 
Garnett, D. R. 1992, AJ, 103, 1330 



Garnett, D. R., Kennicutt, R. C., & BresoUn, P. 2004, ApJ, 607, 
L21 

Gavilan, M., Buell, J. P., MoUa, M. 2005, A&A, 432, 861 
Giavalisco, M. et al. 2004, ApJ, 600, 93 

Gonzalez Delgado, R. M., Leitherer, C., & Heckman, T. 1999, 

ApJS, 125, 489 
Gordon, D. & Gottesman, T. S. 1981, AJ, 86 161 
Groves, B., Dopita, M., & Sutherland, R. 2004, ApJS, 153, 9 
Hillier, D. J., & Miller, D. L. 1998, ApJ, 496, 407 
Hirschi, R., Chiappini, C, Meynet, G., Maeder, A., & Ekstrom, S. 

2008, lAU Symp. 250, Massive Stars as Cosmic Engines, ed. F. 

Bresolin, P. A. Crowtlier, & J. Puis (Cambridge: CUP), 217 
Huehra, J. P., Davis, M., Latham, D., & Tonry, J. 1983, ApJS, 52, 

89 

Hunter, D. A. & Gallagher, J. S., HI. 1986, PASP, 98, 5 

Izotov, Y. I., Stasiriska, G., Meynet, G., Guseva, N. G., & Thuan, 

T. X. 2006, A&A, 448, 955 
Izotov, Y. I., & Thuan, T. X. 2004, ApJ, 616, 768 
Jansen, R. A., Fabricant, D., Franx, M., & Caldwell. N. 2000a, 

ApJS, 126, 331 

Jansen, R. A., Franx, M., Fabricant, D., & Caldwell, N. 2000b, 

ApJS, 126, 271 
Kauffmann, G. et al. 2003, MNRAS, 346, 1055 
Kennicutt, R. C. 1998, ARA&A, 36, 189 

Kewley, L. J., Brown, W. R., Geller, M. J., Kenyon, S. J., & Kurtz, 

M. J. 2007, AJ, 133, 882 
Kewley L. J. & Dopita, M. A. 2002, ApJS, 142, 35 
Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & 

Trevena, J. 2001, ApJ, 556, 121 
Kewley L. J., Ellison, S. L. 2008, ApJ, 681, 1183 
Kewley, L. J., Groves, B., Kauffmann, G., & Heckman, T. 2006, 

MNRAS, 372, 961 
Kewley, L. J., Heisler, C. A., Dopita, M. A., Sutherland, R., Norris, 

R. P., Reynolds, J., & Lumsden, S. 2000, ApJ, 530 704 
Kewley L. J., Jansen, R. A., & Geller, M. K. 2005, PASP, 117, 227 
Kinney, A. L. et al. 1993, ApJS, 86, 5 
Kobulnicky, H. A. & Kewley L. J. 2004, ApJ, 617, 240 
Kong, X. & Cheng, F. Z. 2002, A&A, 389, 845 
Kong, X., Cheng, F. Z., Weiss, A., Chariot, S. 2002, A&A, 396, 

503 

Kudritzki, R. P., Pauldrach, A., Puis, J., & Abbott, D. C. 1989, 

A&A, 219, 205 
Langer, N. 1989, A&A, 220, 135 

Lee, J. C, Kennicutt, R. C, Jose G. Funes, S. J., Sakai, S., & 
Akiyama, S. 2009, ApJ, 692, 1305 



12 



Levesque et al. 



Leithcrcr, C. 2008, lAU Symp. 255, Low-Metallicity Star 

Formation, cd. L. Hunt, S. Madden, & R. Schneider, 305 
Leitherer, C, Schaerer, D., Holdader, J. D., Delgado, R. M.G., 

Robert, C, Kune, D. F., de Mello, D. F., Devest, D., & 

Heckman, T. M., et al. 1999, ApJS, 123, 3 
Lejeune, Th., Cuisinier, F., & Buser, R. 1997, A&AS, 125, 229 
Lejeune, Th., Cuisinier, F., & Buser, R. 1998, A&AS, 130, 65 
Levesque, E. M., Massey, P., Olsen, K. A. G., Plez, B., Meynet, 

G., & Maeder, A., 2006, ApJ, 645, 1102 
Maeder, A. & Conti, P.S. 1994, ARA&A, 32, 227 
Magris, C. G., Binette, L., & Bruzual, A. G. 2003, ApJS, 149, 313 
Mallery, R. P. et al. 2007, ApJS, 173, 482 

Martm-Manjon, M. L., MoUa, M., Di'az, A. I., Terlevich, R. 2008, 

MNRAS, 385, 854 
Matteucci, F. 1986, MNRAS, 221, 911 
Mattcucci, F. & Tosi, M. 1985, MNRAS, 217, 391 
McGaugh, S. 1991, ApJ, 380, 140 

Meynet, G. 1993, in The Feedback of Chemical Evolution on The 

Stellar Content of Galaxies, eds. D. AUoin, G. Stasinska, Publ. 

dc rObservatoirc de Paris, p. 40 
Meynet, G., Maeder, A., Schaller, G., Schaerer, D., & Charbonnel, 

C. 1994, A&AS, 103, 97 
Modjaz, M., Kewley, L. J., Kirshncr, R. P., Stanek, K. Z., Challis, 

P., Garnavich, P. M., Greene, J. E., Kelly, P. L., Prieto, J. L. 

2008, AJ, 135, 1136 
Moustakas, J., Kennicutt, R. C, & Tremoni, C. A. 2006, ApJ, 642, 

775 

Moy, E., Rocca-Volmcrangc, B., & Fioc, M. 2001, A&A, 365, 347 
Noeskc, K. G. et al. 2007a, ApJL, 660, 47 
Noeskc, K. G. et al. 2007b, ApJL, 660, 43 

Osterbrock, D.1989, Astrophysics of gaseous nebulae and a<;tive 

galactic nuclei (University Science Books) 
Pagcl, B. E. J., et al. 1979, MNRAS, 189, 95 

Panuzzo, R., Bressan, A., Granato, G. L., Silva, L., & Danese, L. 

2003, A&A, 409, 99 
Pauldrach, A. W. A., Hoffman, T. L., & Lennon, M. 2001, A&A, 

375, 161 

Peimbert, M. 1967, ApJ, 150, 925 



Pettini, M. & Pagcl, B. E. J. 2004, MNRAS, 348, 59 
Pilyugin, L. S. 2000, A&A, 362, 325 
Ranch, T. 1997, A&A, 320, 237 
Rigby, J. R. & Rieke, G. H. 2004, ApJ, 606, 237 
Rosa-Gonzalez, D., Terlevich, E., & Terlevich, R. 2002, MNRAS, 
332, 283 

Salpeter, E. E. 1955, ApJ, 121, 161 

Schaerer, D., Charbonnel, C, Meynet, G., Maeder, A., & Schaller, 

G. 1993a, A&AS, 102, 339 
Schaerer, D., Meynet, G., Maeder, A., & Schaller, G. 1993b, A&AS, 

98, 523 

Schaerer, D. & Vacca, W .D. 1998, ApJ, 497, 618 
Schaerer, D. & Vacca, W. D. 1998, ApJ, 497, 618 
Schaller, G., Schaerer, D., Meynet, G., Maeder, A. 1992, A&AS, 
96, 269 

Schmutz, W., Leitherer, C, & Gruenwald, R. B. 1992, PASP, 104, 
1164 

Snijders, L., Kewley, L. J., & van der Werf, P. P. 2007, ApJ, 669, 
269 

Stanek, K. Z., et al. 2006, Acta Astron. 56, 333 

Stasinska, G. & Leitherer, C. 1996, ApJS, 107, 661 

Stasinska, G., Fernandes, Cid Fernandes, R., Mateus, A., Sodre, 

L., & Asari, N. V. 2006, MNRAS, 371, 972 
Sutherland, R. S. 1998, MNRAS, 300, 321 
Sutherland, R. S. & Dopita, M. A. 1993, ApJS, 88, 253 
Thuan, T. X. et al. 1999, A&AS, 139, 1 
Thuan, T. X. & Martin, G. E. 1981, ApJ, 247, 923 
Tremonti, C. A. et al. ApJ, 613, 898 

van Zee, L., Salzer, J. J., Haynes, M. P., O'Donoghue, A. A., & 

Balonek, T. J. 2005, AJ, 116, 2805 
Vazquez, G. A. & Leitherer, C. 2005, ApJ, 621, 695 
Vazquez, G. A., Leitherer, C, Schaerer, D., Meynet, G., & Maeder, 

A. 2007, ApJ, 663, 995 
Veilleux, S. & Osterbrock, D. E. 1987, ApJS, 63, 295 
Wirth, G. D. et al. 2004, AJ, 127, 3121 

Yin, S. Y., Liang, Y. C, & Zhang, B. 2007, ASPC, 373, 686 
Zaritsky, D., Kennicutt, R. C, & Huchra, J. P. 1994, ApJ, 420, 87 



Modeling Star-Forming Galaxies 



13 




1000 



\ (Angstrom) ), (Angslrom) 



Fig. 1. — Left: FUV spectra generated by the Starburst99 code adopting an instantaneous burst SFH. The spectra were generated 
using the Geneva HIGH evolutionary models, showing the progression of the spectra with age in 1 Myr increments for the full range of 
metallicities. At 5 Myr the FUV spectra generated by Starburst99 assuming an continuous SFH are also shown (dotted line). Ionization 
potentials for the relevant elements are marked on the x-axis. Right: Subtraction of the FUV spectra when adopting the Geneva HIGH 
evolutionary tracks and the STD evolutionary tracks; HIGH — STD is plotted. In both panels the wavelengths are plotted on a log scale. 
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Fig. 2. — Relative ionization fractions for emission lines produced by Mappings III, plotted as a function of relative distance from inner 
surface of the nebula. The emission features shown here are [Nil] (dashed green line), [Oil] (dashed red line), [OIII] (dashed-dotted red 
line), [SII] (dashed blue line), and [SIII] (dashed-dotted blue line). For these outputs a zero-age instantaneous burst SFH, a metallicity of 
Z = 0.05^0, an ionization parameter q = 1 X 10^, and an electron density rie = 100 cm~^ is assumed. 
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Fig. 3. — Evolution of the diagnostic emission line ratios with age, shown for all five metallicies and ranging from to 5 Myr assuming 
an instantaneous burst star formation history. Models generated with the Geneva HIGH (solid line) and Geneva STD (dashed line) 
evolutionary tracks are compared. An rie = 100 is assumed. 
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Fig. 4. — [NII]/Ha vs. [0III]/H/3 diagnostics for the instantaneous burst SFH model grids evolving from 0.0 Myr to 5.0 Myr in increments 
of 1.0 Myr. The models are plotted with lines of constant metallicity vs. lines of constant ionization parameter. Grids generated with the 
Geneva HIGH evolutionary tracks are plotted with solid lines, while grids generated with the Geneva STD tracks are plotted with dashed 
lines. An electron density rie = 100 is assumed. The grids are compared to our sample of 60,920 SDSS star-forming galaxies from Kewley 
et al. (2006) (points), 95 NFGS galaxies from Jansen et al. (2006b) (blue triangles), blue compact galaxies from Kong & Cheng (2002) (red 
circles), and 10 metal-poor galaxies from Brown et al. (2006) (large open circles). 
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Fig. 5. — [Nil] /Ho vs. [0III]/H/3 diagnostics for the continuous SFH model grids at an age of 5.0 Myr. The models are plotted with lines 
of constant mctallicity vs. lines of constant ionization parameter. Grids generated with the Geneva HIGH evolutionary tracks are plotted 
with solid lines, while grids generated with the Geneva STD tracks are plotted with dashed lines. An electron density = 100 is assumed. 
The grids are compared to our sample of 60,920 SDSS star-forming galaxies from Kewley et al. (2006) (points), 95 NFGS galaxies from 
Jansen et al. (2006b) (blue triangles), blue compact galaxies from Kong & Cheng (2002) (red circles), and 10 metal-poor galaxies from 
Brown et al. (2006) (large open circles). Lines of constant metallicity and ionization parameter follow the legend shown in Figure 4. 
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Fig. 6. — [NII]/[OII] vs. [OIII]/[OII] diagnostics for the instantaneous burst SFH model grids evolving from 0.0 Myr to 5.0 Myr in 
increments of 1.0 Myr. The models are plotted with lines of constant metallicity vs. lines of constant ionization parameter. Grids 
generated with the Geneva HIGH evolutionary tracks are plotted with solid lines, while grids generated with the Geneva STD tracks are 
plotted with dashed lines. An electron density = 100 is assumed. The grids are compared to our sample of 60,920 SDSS star-forming 
galaxies from Kewley et al. (2006) (points), 95 NFGS galaxies from Jansen et al. (2006b) (blue triangles), blue compact galaxies from Kong 
& Cheng (2002) (red circles), and 10 metal-poor galaxies from Brown et al. (2006) (large open circles). Lines of constant metallicity and 
ionization parameter follow the legend shown in Figure 4. 
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Fig. 7. — [NII]/[OII] vs. [OIII]/[OII] diagnostics for the continuous SFH model grids at an age of 5.0 Myr. The models are plotted with 
lines of constant mctallicity vs. lines of constant ionization parameter. Grids generated with the Geneva HIGH evolutionary tracks are 
plotted with solid lines, while grids generated with the Geneva STD tracks are plotted with dashed lines. An electron density rie = 100 is 
assumed. The grids are compared to our sample of 60,920 SDSS star- forming galaxies from Kewley et al. (2006) (points), 95 NFGS galaxies 
from Jansen et al. (2006b) (blue triangles), blue compact galaxies from Kong & Cheng (2002) (red circles), and 10 metal-poor galaxies 
from Brown et al. (2006) (large open circles). Lines of constant metallicity and ionization parameter follow the legend shown in Figure 4. 
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Fig. 8. — [SII]/Ha vs. [OIII]/H/3 diagnostics for the instantaneous burst SFH model grids evolving from 0.0 Myr to 5.0 Myr in increments 
of 1.0 Myr. The models are plotted with lines of constant metallicity vs. lines of constant ionization parameter. Grids generated with the 
Geneva HIGH evolutionary tracks are plotted with solid lines, while grids generated with the Geneva STD tracks are plotted with dashed 
lines. An electron density rie = 100 is assumed. The grids are compared to our sample of 60,920 SDSS star-forming galaxies from Kewley 
et al. (2006) (points), 95 NFGS galaxies from Jansen et al. (2006b) (blue triangles), blue compact galaxies from Kong & Cheng (2002) (red 
circles), and 10 metal-poor galaxies from Brown et al. (2006) (large open circles). Lines of constant metallicity and ionization parameter 
follow the legend shown in Figure 4. 
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Fig. 9. — [SII]/Ho vs. [OIII]/H/3 diagnostics for the continuous SFH model grids at an age of 5.0 Myr. The models are plotted with lines 
of constant metallicity vs. lines of constant ionization parameter. Grids generated with the Geneva HIGH evolutionary tracks are plotted 
with solid lines, while grids generated with the Geneva STD tracks are plotted with dashed lines. An electron density = 100 is assumed. 
The grids are compared to our sample of 60,920 SDSS star-forming galaxies from Kewley et al. (2006) (points), 95 NFGS galaxies from 
Jansen et al. (2006b) (blue triangles), blue compact galaxies from Kong & Cheng (2002) (red circles), and 10 metal-poor galaxies from 
Brown et al. (2006) (large open circles). Lines of constant metallicity and ionization parameter follow the legend shown in Figure 4. 
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Fig. 10. — Left: Comparison of instantaneous burst model grids from Kewley et al. (2001; green dash-dotted lines), Dopita et al. (2006) 
(blue; dashed hnes), and this work (red sohd lines) for the [NII]/Ha vs. [OIII]/H/3 (top), [NII]/[OII] vs. [OIII]/[OII] (middle), and [SII]/Ha 
vs. [0III]/H/3 (bottom) diagnostic diagrams. The Kewley et al. (2001) diagnostics range from g = 1 X 10^ to g = 3 X 10*, and use the 
Starburst99 stellar population synthesis models and the Geneva evolutionary tracks. The Dopita et al. (2006) grids range from 0.0 to 4.0 
Myr in increments of 1.0 Myr. The models are plotted with the SDSS galaxies (Kewley et al. 2006) and metal-poor galaxies (Brown et al. 
2006). Right: Comparing the 8.0 Myr Starburst99/Geneva model grids from Kewley et al. (2001) to the 5.0 Myr continuous SFH models 
from this work (red solid lines). The models are plotted with the SDSS galaxies (Kewley et al. 2006) and metal-poor galaxies (Brown et 
al. 2006). 
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Fig. 11. — Comparison of the FUV ionizing spectrum generating by Starburst99 when adopting the Geneva HIGH tracks (solid black 
line), the Geneva STD tracks (dashed red line), and the newest generation of the Geneva tracks which include the effects of rotation 
(dashed-dotted blue line). It is apparent that the rotating tracks generate a much harder ionizing spectrum, particularly in the high-energy 
regime. All of these tracks are at a metallicity o{ Z = Zq. An age of 5.0Myr, a continuous SFH, and rie = 100 are assumed. 



